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I. INTRODUCTION 

Understanding the dynamics of quantum fields in 
curved spacetime is a topic of general interest with im¬ 
portant applications to early Universe cosmology or black 
hole physics. Prominent examples are the primordial 
density fluctuations during inflation [ij and the Hawk- 
ing/Unruh radiation from (analog) black holes M. The 
case of de Sitter space has attracted particular attention 
both because of its relevance for inflationary cosmology 
and because it provides a simple, maximally symmetric 
example where the nontrivial aspects of the curved ge¬ 
ometry come into play. 

The investigation of quantum field dynamics on de Sit¬ 
ter space has many facets. A general class of studies con¬ 
cerns the solutions of the field equations of motion in a 
fully de Sitter symmetric state. Renormalizability and 
the equivalence principle select the so-called Chernikov- 
Tagirov-Bunch-Davies (CTBD) state [^, @ out of the 
class of Allen-Mottola a-vacua Topical questions 

include the calculation of radiative corrections to the field 
dynamics [Tol - E^ , the relevance of global coordinate sys¬ 
tems as compared to the (expanding) Poincare patch, of 
interest for inflationary cosmology [29l - [^ . and the rela¬ 
tion between Lorentzian and Euclidean de Sitter spaces 
. Other studies concern other (possibly non de Sit¬ 
ter invariant sta tes, e.g., in the context of inflationary 
cosmology 0, [13, , or in relation with the question of 

the quantum stability of de Sitter space 0,[13, |3l|,[H, [H- 

lil- 

An important issue is the understanding of interacting 
field theories. Not only are these more difficult to tackle 
than in flat space for obvious technical reasons, but the 
nontrivial gravitational field leads to specific effects with 
no flat space analog. Striking examples are the possibility 


of curvature-induced phase transitio ns [R)I-I1^ , the decay 
of massive particles into themselves 1^ , th e genera¬ 
tion of a nonvanishing photon mass [l3l l45l l46l| , and the 
radiative restoration of spontaneously broken symmetries 
[13, . In the CTBD vacuum, light scalar fields in 

units of the spacetime curvature exhibit large, gravita¬ 
tionally enhanced quantum fluctuations on superhorizon 
scales. Notably, this is at the origin of the scale in¬ 
variance of the (tree-level) power spectrum of primordial 
density fluctuations in inflationary cosmology. However, 
such large fluctuations are also responsible for the strong 
infrared sensitivity of radiative corrections. Perturba¬ 
tion theory typically exhibits infrared and secular (large 
time) divergences flsj - flTi The former are generic 
for massless bosonic field theories [13] while the latter 
are characteristic of nonequilibrium^ systems [53| . Both 
types of (spurious) divergences signal a breakdown of the 
perturbative expansion and need to be resummed. 


Various resummation methods or genuine nonpertur¬ 
bative approaches have been devised to deal with such 
issues in de Sitter space. The most prominent one is the 
stochastic effective theory put forward in Refs. [s^. [53|. 
This has been shown to correctly capture the nonper¬ 
turbative dynamics of superhorizon modes at lea ding 
order accuracy both in secular logarithms OSS Hi 
and in infrared enhancement factors [13, US]. More re¬ 
cent approaches employ various quantum field theoreti¬ 
cal tools such as large-V techniques 
renormalization group methods [isTiinMnej the 


^ The nonequilibrium nature of (the Poincare patch of) de Sitter 
space stems for the cosmological expansion in standard comoving 
coordinates. 
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Wigner-Weisskopf approach [H, , reduced density ma¬ 

trices [g^, and solutions of the Kadanoff-Baym equa¬ 
tions (KBEs) [g^ Ig^ — the nonequilibrium version of the 
D^on-Sch wi^ey quations—and their Boltzmann limit 

In Ref. [Z^, we have studied the KBEs for the cor¬ 
relator of an 0{N) scalar field theory with quartic self¬ 
interactions using the physical momentum representation 
of de Sitter correlators [76l - [7§| . We have obtained an ex¬ 
act analytical solution for the two-point correlator in the 
regime of superhorizon momenta when the self-energy is 
computed at two-loop order in perturbation theory. The 
secular di verg ences of perturbative calculations of the 
correlator |15l - ll7l l28l| are a mere artifact of expanding 
the latter in terms of self-energy insertions at a finite or¬ 
der. Solving the KBEs amounts to resumming the infinite 
series of such self-energy insertions—and thus the associ¬ 
ated large infrared logarithms—which results in a modi¬ 
fied power law behavior of the correlator. This is similar 
to the generation of an anomalous dimension for critical 
systems in statistical physics. From the (resumed) prop¬ 
agator, we could compute various quantities of interest, 
such as the infrared field strength and mass renormaliza¬ 
tion, at two-loop order. The applicability of such pertur¬ 
bative computations is limited to not too light fields. The 
cases of fields with either vanishing or negative tree-level 
square mass suffer from infrared divergences and require 
further resummations. 

In the present article, we generalize this approach to 
the nonperturbative 1/N expansion at next-to-leading 
order (NLO). The leading-order (LO) approximation re¬ 
sums the infrared divergences of perturbation theory into 
a self-consistent (local) mass term. This captures inter¬ 
esting nontrivial physics, such as dynamical mass gen¬ 
eration and radiative symmetry restoration |4^, or 
nonperturbative quantum contributions to non-Gaussian 
correlators [2^, [g^. The NLO approximation involves 
an infinite series of nonlocal multiloop contributions to 
the self-energy. For infrared momenta, the latter can be 
summed in a closed form using the results of Ref. [g^. 
The NLO self-energy has a similar form as the two-loop 
expression, however, with different infrared power law 
exponents. This allows us to obtain an exact analyti¬ 
cal solution of the corresponding KBEs in the infrared 
regime using the techniques developed in our previous 
work 

We thus obtain the complete spacetime structure of the 
propagator for superhorizon momenta at NLO in the 1 /N 
expansion. As in the two-loop case, the infrared loga¬ 
rithms of perturbation theory resum into a superposition 
of well-behaved momentum power laws, which charac¬ 
terize the decay of the field correlator at large spacetime 
separations. We compute the leading infrared behavior 
for deep superhorizon momenta, from which we extract 
the field strength and mass renormalization. These re¬ 
ceive NLO corrections which are nonperturbative func¬ 
tions of the field self-coupling. We study the cases of van¬ 
ishing and negative tree-level square mass, which both 


correspond to effectively strongly coupled regimes in the 
infrared. Finally, we compute the local field variance, 
which agrees with the result of the stochastic approach 
at the same approximation order. 

The paper is organized as follows. Section HIl describes 
the setup and briefly reviews the formulation of KBEs 
in the p-representation [t^ . We compute the self-energy 
at NLO in the 1/N expansion in Sec. Illll and present the 
solution of the corresponding KBEs in Sec. lIVI The NLO 
result requires the solution of an appropriate gap equa¬ 
tion for the local mass, which is described in Sec. El We 
discuss our analytical solution for the field correlator in 
various regimes of interest in Sec. ED and we conclude in 
Sec. IVIII Finally, some technical details and calculations 
are presented in Appendixes lADCl 

II. GENERAL SETTING: KBES IN THE 
p-REPRESENTATION 

We consider an 0{N) scalar field theory on the ex¬ 
panding Poincare patch of de Sitter space in D = d + 1 
dimensions. In conformal time — oo < ry < 0 and co¬ 
moving spatial coordinates X, the invariant line element 
reads (we set the Hubble scale H = 1) 

ds^ = r)-^{-drf + dX.dX). (1) 

The classical action is given by 

S = I (□ - TOds) ‘fa - -^(.fafaf^ , (2) 

with = f d^x^—g{x) the invariant measure— g{x) 
is the determinant of the metric—and where summa¬ 
tion over repeated indices a = is understood. 

Here, □ is the appropriate Laplace-Beltrami operator 
and m^g = m? + //R. includes a possible coupling to the 
Ricci scalar TZ = d{d+l). We consider a symmetric state 
such that {ipa) = 0 and the correlator G and the self¬ 
energy S are diagonal, e.g.. Gab = SatG. In the rest of 
the paper we assume a de Sitter invariant state given by 
the free field CTBD vacuum in the remote past ry —>■ — oo. 

The covariant inverse propagator is given by 

G~^{x,x') = Gq^ {x,x') —'E,{x,x'), (3) 

where 

iG//^{x,x')= {na:-mls)S^^\x,x'), (4) 

with 5^^\x,x') = 5^^\x — x')/y/—g{x). Extracting a 
possible local part from the self-energy,^ one writes 

T,{x, x') = {x, x') + E(a:, x')^ (5) 


^ One may have to include more complicated structures in the 
local contribution to the self-energy when discussing ultraviolet 
renormalization. For instance in Z) = 4, there ap pea rs a term 
(fc, fc') due to field-strength renormalization [l^ . 
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where the local a part is constant for a de Sitter invariant 
state. We include it in a redefinition of the mass 

= mls + cr ( 6 ) 

and define, accordingly, the propagator 

iGl^{x,x') = (□, - m2) S^^\x,x'), (7) 

in terms of which we have 

G~^{x,x') = G]^{x,x') —'E{x,x'). (8) 

The aim of the present work is to solve this equation for 
G when the self-energy ([S]) is computed at NLO in a 1/A^ 
expansion. The main difficulty concerns the convolution 
with the nonlocal kernel 51(a;, x') when inverting Eq. ([5|). 
To this aim, we shall exploit the techniques developed in 
Refs. [60, which rely on exploiting the de Sitter 

symmetries in a physical momentum representation, the 
so-called p-representation. 

Exploiting the spatial homogeneity and isotropy in co¬ 
moving coordinates, one writes 

G{x,x')=J (9) 

De Sitter symmetries guarantee that the correlator admit 
the following scaling form [t^ 

( 10 ) 

where p = —Krj and p' = —Krj' are the physi¬ 
cal momenta associated with the comoving momentum 
K at times p and p' respectively. Similarly, the p- 
representation of the self-energy is 

t{p,p',K) = {pp')^K^t{p,p'). (11) 

Solving the Schwinger-Dyson equation ([5]) for the prop¬ 
agator G for a given self-energy D in de Sitter space can 
be viewed as an initial value problem with initial data 
to be specified in the infinite past p —>■ —oo. This can 
be formulated by introducing a closed contour in time— 
the so-called in-in or Schwinger-Keldysh formalism — 

which allows one to conveniently grab together the vari¬ 
ous components of Green’s functions. Exploiting the way 
time and momentum are tight together by gravitational 
redshift in de Sitter spacetime, the time evolution can be 
traded for a momentum evolution, with initial data to be 
specified at p —>■ -l-oo. Introducing a closed contour C in 
momentum, the propagator reads [t^ 

G(p,p') = E(p,p') - ^signc-(p-p')p(p,p') (12) 


P{PtP') = F{p':P)^ and p(p,p') = —p(p',p). The (nonlo¬ 
cal) self-energy S(p,p') admits a similar decomposition. 

The Dyson-Schwinger equations expressed in the in-in 
formalism are called the Kadanoff-Baym equations [6^, 
113 • In the p-representation, they read [t^ 


1 / 2-1 
2,1 4 




poo 

F{p,p')+ ds'Ep{p,s)F{s,p') 

Jp 

poo 

= / dstp{p,s)p{s,p') (13) 


and 


dl + l- 




Kp,p')= dstp{j),s)p{s,p'), 

Jp' 

(14) 


where we introduced (the last equality defines e) 



In the following, we consider light fields in units of the 
spacetime curvature, i.e., e ~ M"^ jd <?; 1. 

As announced, momentum here plays the role of the 
time-evolution variable. It is remarkable that the KBEs 
effectively reduce to a (0 -I- l)-dimensional dynamical 
problem in the p-representation. In fact, the space di¬ 
mensionality is completely hidden in the expression of 
the self-energy through d-dimensional loop integrals. The 
statistical function F encodes the information about 
the actual quantum state of the system. Having in 
mind an adiabatic switch-on of the interactions, the ini¬ 
tial data corresponding to the CTBD vacuum are given 
by F[p,p )|p=p/_>_|_oo = 1/2, dpFijp^p )|p=p/_>_|_Qo = 0, 
i9pi9p'F(p,p')|p=p/^+oo = 1/2. The nontrivial initial data 
for the spectral function is determined by the equal-time 
commutation relations: clp/5(p,p')|p=p' = —1. 


III. THE SELF-ENERGY AT NLO 


We now come to the calculation of the self-energy at 
NLO in a 1/A expansion. This has been discussed in the 
context of nonequilibrium QFT in Minkowski spacetime 
in Refs. IziS- In de Sitter spacetime the corresponding 
expressions have been derived in the p-representation in 
Ref. [t^ . Here, we briefly review the material relevant 
for our present purposes. 

The local contribution in Eq. m can be written in 
closed form in terms of the exact propagator as 



Fiq,q) 

q 


(16) 


where F and p denote the p-representations of the statis¬ 
tical and spectral two-point functions respectively. Here, 
the sign function is to be understood along the contour C; 
see Ref. [t^ for details. Notice the symmetry properties 


where /q = / d‘^g/(27r)'^. This is represented diagram- 
matically in Fig. [TJ The corresponding expression at 
NLO is obtained by plugging the (yet unknown) NLO 
correlator on the right-hand side and discarding terms 
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FIG. 1: The local contribution to the self-energy E(x, x')\ see 
Eq. ©. The internal line in the diagram is given by the full 
propagator G{x,x). 


of 0(1/7V^). As we shall see below, this actually pro¬ 
vides a self-consistent equation for the so-called tadpole 
mass M defined in Eq. ®. For now, we do not need the 
expression of M and leave it as a free parameter in the 
calculation. 



FIG. 2: A typical multiloop diagram contributing to the non¬ 
local self-energy T,{x,x') at NLO. The latter actually resums 
all diagrams of similar topology with an arbitrary number 
n > 1 of bubbles in the upper part, as described by the inte¬ 
gral equation mi) in the p-representation. 


The nonlocal part of the self-energy reads, at NLO, 

S(p,p') = {pp')^ [ -Gm (qp, qp') I (rp, rp') 

3iV q 

(17) 

where r = |e -|- q|, with e an arbitrary unit vector, and 
where the function / resums the infinite series of bubble 
diagrams represented in Fig. [2] It solves the following 
integral equation^: 

i[PiP') = n(P>p') -i J^dsfl{p,s)i{s,p'), (18) 


because the LO propagator is given by Gmo with Mq the 
LO tadpole mass. Here, we shall keep Gm in the inter¬ 
mediate steps of the calculation for notational simplicity 
and consistently expand the final expressions in 1/7V. 

Let us finally mention that H, /, and S all have similar 
decompositions as in Eq. © on the momentum contour 
C. The product GmGm in the expression (HU) of the 
function H gives rise to the combinations FmFm — jPmPm 
for Hp and 2FmPm for Hp. Similarly, the product GmI 
in dm) yields FmIf - jPAiIp for and FmIp + PmIf 
for Sp. The explicit integral equations satisfied by the 
components If and Ip, obtained from Eq. dm, read 

Ipip.p') =%(j>,p') + [ dsflp(j),s)ip{s,p') ( 20 ) 

Jp 

and 

poo 

If{p,p') = fiFip,p') - / dsflF{p,s)ip{s,p') 

Jp' 

pOO 

+ / dsflp{p,s)iF{s,p'). (21) 

Jp 

Let us recall the approximation strategy used in 
Refs. [63, [73 for the calculation of correlators for su¬ 
perhorizon momenta ^ 1. For light fields in units of 
the spacetime curvature, the relevant dynamics is domi¬ 
nated by the gravitationally amplified superhorizon fluc¬ 
tuations, and we shall neglect the effect of interactions 
for subhorizon modes. We thus consider a free evolu¬ 
tion up to a given scale ^ 1 below which we fully 
take into account the field self-interaction. In practice, 
this amounts to restricting the convolution integrals^ in 
Eqs. ini) and dUl) to and We ob¬ 

tain a closed set of integro-differential equations which 
only involve two-point functions for superhorizon mo¬ 
menta. One can show, using the method developed in 
Ref. [^, that the dominant contribution to the loop 
integrals in Eqs. dm and dm for p,p' < pL comes from 
loop momenta qp,qp',rp,rp' < p. These integrals can 
thus be evaluated by using the leading infrared behavior 
of the propagator, that is. 


where the momentum convolution on the right-hand side 
is to be taken on the momentum contour C to account 
for the in-in formulation of the problem and where H is 
the one-loop bubble integral 

fV/ /X A. f Gm {qp,qp') Gm {rp,rp') 

n(p,p =--(pp 2 /-^ 

Q Jq q r 

(19) 

with r = |e-|-q|. Notice that we have used the propagator 
Gm in the NLO expressions dm-®- This is consistent 


Fm{p,p') = \/^j^, (22) 

P^m{p,p') = -\fpp'Fi, ^In , (23) 

where F^, = /47r and Vviu) = sivl\i{vu) / v . The 

one-loop bubble integral (fTOl) has been computed in 
and the solution of the integral equation Eq. (TT^ for in¬ 
frared momenta has been obtained in closed form. The 


^ This is directly related to the Dyson-Schwinger equation for the 
four-point vertex in the large-limit [goII . 


^ The role of subhorizon modes in Eqs. 1201 and has been 

studied in detail in Ref. [ 6 Q|| . They modify the result m and 
112511 below by a simple renormalization factor of order unity. 
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dominant infrared behaviors of the statistical and spec¬ 
tral components of the infinite series of bubbles I read 


i¥^ip,p') 




Vi¥ {pp'Y 


VpY 


n 


->-e ’ 

(24) 

P' ) ’ 

(25) 


where e has been defined in Eq. (USD, 


_ fid A 

^ ~ 6e {2 t:Y SM^flo+i ’ 


(26) 


with fid = 27r'^/^/r(c?/2), id = ^/lF — TTp, and where we 
defined 'P^{u) = Va{u)e~^'^'^'^. Equations (IMl) and (1^ 
are valid provided TTp <g; 1. More general expressions can 
be found in [i^- Writing v = d/2 — e, we thus have 

£«£-!- TTp/d < 1. 

As discussed in Ref. [^, the function / is of the very 
same form as the one-loop bubble 11 with the replacement 
Id ^ D. Formally expanding Eqs. (El and (El) in tt p 
generates the infinite series of bubble diagrams, each of 
which brings an additional power of infrared (secular) 
logarithms TTp In(pp') and 7 rpln(p/p'). The infinite series 
of such infrared logarithms resum into the modified power 
laws (El and (1^ . 

Using the expressions (0^ - (l^l) . we can evaluate the 
dominant infrared behavior of the loop integrals in (1171) . 
This is detailed in Appendix lAl We obtain 


^f{p,p') 


Op Fjj 

(pp')3/2 ’ 


•p27 

(pp')3/2 ^ 



(27) 

(28) 


where we defined 7 = (e -I- £)/2 and 


<j 


p ~ 


p 

N 



(29) 


We emphasize that, if it is justified to neglect rela¬ 
tive corrections of order £, £, and 7 in numerical pref¬ 
actors, one should not do so in the exponents of the 
various power law dependences in the momenta. In¬ 
deed, such corrections become relevant at large values of 
I ln(ip/ip')| involved in convolution integrals, as discussed 
in Refs. [63,[l3|; see also 


IV. SOLUTION OF THE KBES IN THE 
INFRARED 

It is remarkable that the NLO self-energy, Eqs. El 
and (El! has the very same form® as the two-loop one 


® More precisely, we obtain the same expressions as Eqs. (19) and 
(20) in Ref. with s(u) = and cf{u) = V/P{u). We 


obtained in [74l | up to the expression of Op and with the 
replacement £ —>■ 7 . This is rooted in the fact that the 
two-loop self-energy can be written as oc oc 

IIGm, with n the one-loop bubble defined in Eq. m- 
Now, the NLO self-energy assumes a similar form with 
the one-loop bubble replaced by the infinite series of bub¬ 
ble diagrams, that is, 11 —>■ /; see Eq. El- As men¬ 
tioned above this simply amounts to a modified exponent 
V ^ id. The calculation of the loop integral in Eq. El 
is then essentially the same as in the two-loop case; see 
Appendix 

This observation actually allows us to directly use the 
results of our previous work [tJ for the solution of the 
KBEs (TT^ and (11^1) . There, we had shown that, for su¬ 
perhorizon momenta p,p' < p, these reduce to a single 
integro-differential equation for a function of one vari¬ 
able, which can be solved exactly by analytical means. 
Here, we simply quote the relevant results and refer the 
reader to Ref. for details. 

The statistical and spectral components of the corre¬ 
lator read, for p,p' < p, 


F^^{P,P') = VpYF^ 
and 


c+ 


(j)p'Y+- 


1 

, (30) 

c- 

1 

1 

s 

p') ’ 

(31) 


with 


p{u) = c+VY (.u) + C-V2_ {u) 

—7)^) T d——. y(^); 


(32) 


where the approximate expression in the second line is 
valid for |u| > 1. Here, c+ -I- c_ = 1, with, up to correc¬ 
tions® 0{N-Y, 


161/272 4 (i ^72 ’ 


and 


id± = id±± 


8idid±'y 


id± ± 


2d?^^ 


(33) 


(34) 


where id± = z/ ± 7 , and where we neglected terms 
0{e/N,^/N) in the final expressions. 

As we already observed in the two-loop approximation 
[tJ, we find that the NLO propagator (l3^ - (l32|) is es¬ 
sentially described by a linear superposition of two free 
massive field propagators with masses m± defined as 


Id? „ d 

- 7 = V ^ = 2 - e±- 


(35) 


check that by setting 7 —>■ e in the present Eq. m , we recover 
the large-AC limit of the two-loop expression of o-p] see Eq. (21) 
of Ref. [ 3 . 

The actual solution of the KBEs [3 yields c± = {i> ^ za)/(2i>) 
and u± = :k 2'yv -\- 7 ^, where u = -\- apKA.'y'^). 
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Notice that m+ < m_. Expanding the above solution in 
powers of Up formally corresponds to the infinite series 
of nonlocal self-energy insertions and generates arbitrary 
powers of infrared logarithms (Tpln(pp') or tTpln(p/p'). 
These correspond to the infrared secular logarithms ob¬ 
tained in strict perturbative calculations of the field cor¬ 
relator |15l - [r^ . Again, the summation of these con¬ 

tributions through the KBEs yields well-behaved infrared 
power laws with modified exponents v ^ v± — Going 
back to spacetime variables, the NLO correlator can be 
written 


uated at LO. Equation (p 8 )) is solved as 

2 Aeff (l + Aeff + Agjj) 
(1 + Aeff)2 

where 


(39) 






612 


D+l 


(40) 


is the LO solution (notice that Mq > 0) and where 
we introduced the effective coupling 


G{x,x') = C+Grn+iz) +C-G,ri-{z), (36) 

where Gm(z) is the propagator of a free scalar field with 
mass^ m <C 1 and z = z(x,x') is the de Sitter invariant 
distance between the points x and x'; see Appendix [BJ 
The solution dsni)-® is expressed in terms of the 
tadpole mass M. The latter has to be determined self- 
consistently from the gap equation ©, with the expres¬ 
sion (fTBl) . 


AefF pp, 11^4 • ( 41 ) 

The relevance of the latter stems from the observa¬ 
tion that the effective potential for light scalar fields 
on superhorizon scales is described by an effective zero¬ 
dimensional field theory [s^, [ll|. Equation (HT|) corre¬ 
sponds to the dimensionless effective coupling of such a 
theory at LO in the 1 /N expansion. 


V. DETERMINATION OF THE TADPOLE 
MASS AT NLO 

The tadpole diagram (1161) involves the correlator in 
the coincidence limit® G{x,x) = {ip^{x))/N. The latter 
is dominated by the strongly amplified superhorizon fluc¬ 
tuations and can be evaluated using the expression (1501) . 
This yields 

^ f F{p,p) 

N J {2TTy p 

^ Fy^d f c_ A 

{2'Ky \d — 2 (P+ — 7 ) d — 2 (f'_ — 7 )/ 

~- - - ( 1 H- — -bO . 

no+iM^ \ N 2d^ej ^ > j 

(37) 


VI. DISCUSSION 

Let us now discuss some of the peculiar features of 
the NLO propagator obtained above. In particular, we 
consider the regime of deep infrared momenta and the 
corresponding behavior of the field correlator at large 
spacetime separations as well as the field variance, given 
by the correlator in the coincident limit. The latter can 
be compared to the results from other methods such as 
the stochastic approach or the Euclidean de Sitter ap¬ 
proach. We shall examine the relevant quantities in var¬ 
ious regimes of parameters from the case of light but 
massive fields, where perturbation theory applies, to the 
cases of either zero or negative tree level square mass, 
described by an effective strongly coupled regime. 

A. Deep infrared momenta and large spacetime 
separations 


The NLO gap equation (l 6 |) thus reads 




™dS + 


A 

6 f2D+iM2 


1 


2 L , ( 70 -^ 0 )^ 
N V £o7o 


(38) 


where we have used Tip/ [2d) ~ 7 — £ and where the sub¬ 
script 0 indicates that the corresponding quantity is eval- 


^ Recall that the present definition of the square mass includes a 
possible nonminimal coupling to the Ricci scalar. 

® The ultraviolet divergence of the latter can be absorbed by a 
standard renormalization of the parameters of the Lagrangian 
m?, and A in D < 4 . The choice of CTBD vacuum guarantees 
that the required counterterms are the same as in Minkowski 
space; see, e.g., 1221 l6Cll. I8III . We shall disregard these aspects 
here and focus on the IR contribution. 


The dominant contribution to the regime of deep in¬ 
frared momenta <C p and/or large momentum (time) 
separation | \n[p/p')\ oc — ^ 1 [where t = — ln(— 77 ) is 

the cosmological time] is governed by the term with the 
lowest mass m+ on the right-hand side of either Eq. (1501) 
or (1551) . In this regime, the NLO propagator reduces to 
that of a renormalized massive field 

(42) 

~-\/pp'^r’.iR (in^l , (43) 

where Z = c+ = 1 — c- and, denoting m\ = 


VlK = 



(44) 
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The field renormalization factor is completely due to 
the NLO nonlocal self-energy and reads, in terms of the 
effective coupling (IdTll . 


Z = 1 - 




1 - 


1 2Aeff) 

^ (1 + Aeff)3 ■ 


(45) 


Similarly, the contribution to the infrared mass due to 
nonlocal self-energy insertions is 


m 


IR 


M2 


= 1 - 


2d‘^e-f 


1 - 


1 A2g(3 + 2Aeff) 
N (1 + Aeff)2 ■ 


(46) 



Taking into account the NLO tadpole contribution, we 
finally obtain 


™IR _ 1 I 1 "^eff (2 Aeff) , , 

M2 (1 + Aeff)2 • 

We see that the factor Z < 1 controls the amplitude 
of the deep infrared fluctuations in momentum space, 
whereas rrij^ controls the momentum dependence of the 
correlator and, in turn, its behavior at large spacetime 
separation, as described in Appendix |B] 


B. Field variance 


Another quantity of interest is the field correlator in 
the coincidence limit, which measures the local field vari¬ 
ance, G{x,x) = {(p‘^{x))/N. The latter can be character¬ 
ized in terms of an effective mass mdyn, defined as (^.1^ 


N ~ riD+iml^J 


(48) 


in analogy with the expression for a free massive field. 
In fact, such a square mass is nothing but the curvature 
of the effective potential of the dimensionally reduced 
theory [5l| . We have already evaluated the field variance 
in Eq. (Eg. This yields, for the contribution of nonlocal 
self-energy insertions. 


^lyn , 2 A2ff 

M2 Nl + X^g' 


(49) 


Including the contribution from the NLO tadpole mass, 
we get 


^lyn ^ 2 Aeff 

M2 +N(l + Aeff)2- 


(50) 


This expression can be compared to the result of the 
stochastic approach, where the local field fluctuations 
can be computed at all orders of perturbation theory 
with leading infrared logarithmic accuracy [i3- For light 
fields, this is equivalently given by the variance of the 
zero mode on Euclidean de Sitter space [1^ . We present 
the calculation of the field variance at NLO in the 1/A^ 


FIG. 3: The dimensionless quantities (from bottom to top) 
Z, rrii^/MQ, rrigy^lMQ, and mjj^/(ZMQ) as functions of the 
dimensionless effective coupling Aeff for N = 4. The two 
lower curves asymptote to 1 — 1/N, whereas the two upper 
ones asymptote to 1. Deviations from unity are due to both 
local and nonlocal 1/A corrections. 


expansion in the stochastic and Euclidean de Sitter ap¬ 
proaches in Appendix [Cl which completely agree with 
Eq. (I50l) . 

We emphasize that the result dSOl) differs from what 
one would obtain by evaluating the field variance only 
from the deep IR behavior, Eq. ([55]). Indeed, this would 
lead to {(p‘^{x))/N —>• The term oc c_ in 

Eq. (l37ll is needed to get the correct variance. However, 
we mention that the ratio turns out to be numer¬ 

ically close to TO^yn arbitrary Aeff. The dimensionless 
quantities Z, mjj^/Mg, ™ir/(^Mo2 ), and m\JMl are 
plotted against the effective coupling AeS in Fig. [3] No¬ 
tice the hierarchy < wliy^IZ. It is important 

to notice that all these square masses are equal to Mg 
at LO and that their differences are entirely due to the 
nonlocal self-energy. Since neither the LO nor the NLO 
tadpole square masses Mg and M'^ are observable, it is 
interesting to consider ratios of physical quantities such 
as 


m 


■IR 


Oyn 


m 


IR 




. 1 A2ff 

N{1 + Aeff)2 ’ 

1 + A -- 

^ 2A(1 + Aeff)3’ 


(51) 

(52) 


whose deviation from unity is entirely an effect of re¬ 
suming the large infrared logarithms from nonlocal self¬ 
energy insertions. The ratios (l5T|) and (1521) are plotted 
against Aeff in Fig. [H 

Finally, we note that all the NLO corrections computed 
in this sections, Eqs. (1451) . (l47ll . (15(1 . (1511 . and (l52ll . ex¬ 
hibit a nontrivial dependence in the (effective) coupling 
and are actually bounded functions of the latter due to 
nonperturbative 1/A corrections. The NLO results thus 
remain valid for arbitrarily large effective couplings. We 
shall now discuss the various cases of interest, from the 
























FIG. 4: The dimensionless ratios of physical quantities 
^m/^dyn (lower curve) and (upper curve) 

as functions of the dimensionless effective coupling Aefi for 
N = 4. These asymptote to 1 — 1/N and 1, respectively 
at large coupling. Deviations from unity are entirely due to 
nonlocal 1/A^ corrections. 





FIG. 5: The dimensionless effective coupling eu as a func¬ 
tion of the tree-level square mass m|g for A = 10“^ in 
D — 3 -h 1 dimen sions. Large positive values of m^g in units 
of y/\/(Q^ d+i) ~ 0.025 correspond to perturbatively small 
effective couplings, whereas small or large negative values cor¬ 
respond to Aefi ~ 1 and Aeff 1, respectively. 


weak coupling, perturbative regime to the strongly cou¬ 
pled one. 


C. Perturbative regime 


Perturbation theory makes sense whenever the 
quadratic part of the action gives the dominant con¬ 
tribution to physical observables. This corresponds to 
the case of light but massive field, with y/\ ^ m^g <C 1 
[13, m, [6l|. In that case, one has Mg = mgg[l -h Aq — 
Ag -h C’(Ag)], Aeff = Ao[l - 2 A 0 -h C>(Ag)] < 1, where we 
defined 


(Wds)^ 


We get 


Z = 


1 - 


2N 


(53) 


(54) 


and 


"ilR — 




‘dyn 


MS 


1 + 



(1 + ]^)ao + C’(Ao) , 

(55) 

(1 + +^(Ao) • 

(56) 


These expressions agree with the two-loop results of 
Ref. [t^ at the relevant order in 1/N, as they should.® 


D. Massless case 


More interesting for our present purposes are theories 
with very small tree-level mass, such that Aq ^ 1, or with 
negative tree-level mass, m^g < 0. Both cases correspond 
to strongly coupled effective theories in the infrared, as 
measured by the effective coupling 6ID; see Fig. [SJ We 
first consider the light mass case. It is well known that, 
even for massless fields, the scalar self-interactions gen¬ 
erate a nonvanishing effective mass. At LO the latter is 
Mg = ^/X/{6nD+l) and the effective coupling AeS ~ 1. 
We thus get, at NLO, 


and 


Z = 1 - 


5 


(57) 


™IR — 

"llyn = 


(' + iv) ■ 
\/«!^ (' + 2v)' 


E. Negative square mass case 

Finally, in the case of a negative tree-level square 
mass, the apparently broken symmetry gets radiatively 
restored by the strong, grav itationally amplified super¬ 
horizon fluctuations |43, |4^. For Ag <C 1, one has 
Mg Ri A/(6fl£)+i|mgg|) = Aglruggl and the effective in¬ 
frared theory is strongly coupled: Aeff ~ 1/Ag 1. We 
obtain 



® The effective coupling introduced in Ref. [tJ] is A = Ao(l-|-2/Af). 


(60) 
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and 


WiR = Aolwdsl (^1-: (61) 

”^dyn = ^ol^T^dsI + ]y^0^ • (62) 

VII. SUMMARY AND PERSPECTIVES 

Quantum field theoretical calculations in curved space- 
times are technically involved as compared to their 
Minkowski counterparts because of the nontrivial sym¬ 
metries of the background geometry. Even the calcula¬ 
tion of the basic two-point correlator of a self-interacting 
scalar field in the maximally symmetric de Sitter space 
is a notoriously difficult task due to the nontrivial con¬ 
volution/memory integrals involving the nonlocal self¬ 
energy. For cosmological spacetimes, the problem can 
be cast in the form of a nonequilibrium, initial value 
setup [S^- A strict perturbative expansion in powers 
of self-energy insertions leads to spurious secular diver¬ 
gences typical of nonequilibrium systems [s^. Solving 
the integro-differential KBEs is formally equivalent to re¬ 
summing the infinite series of self-energy insertions and 
yields a well-defined late-time behavior. Together with 
Refs. [^ . the present work provides techniques to formu¬ 
late and solve KBEs in de Sitter spacetime in the case of 
light fields. Related work for massive fields or for non de 
Sitter invariant states can be found in Refs [H, IH, IH- 

Izi- 

In the present article, we have computed the two- 
point correlator of an 0{N) scalar field in the regime 
of superhorizon momenta at NLO in a 1/N expansion. 
Remarkably the solution can be obtained in closed an¬ 
alytical form. This resums the spurious infrared and 
secular divergences of perturbation theory into a well- 
defined expression. The nonperturbative 1/N expansion 
allows us to compute various quantities of interest, such 
as the infrared field strength and mass renormalization, 
or the local field variance, for arbitrary values of the field 
self-coupling. In particular, this is a valid approxima¬ 
tion scheme for the cases of either vanishing or negative 
tree-level square mass, which correspond to effectively 
strongly coupled regimes not accessible by perturbative 
means. 

The local field variance {(p‘^{x)) obtained here fully 
agrees with the result of either the stochastic, or the Eu¬ 
clidean de Sitter approaches—which are in fact equiva¬ 
lent for what concerns this observable [12,1121 — as well 
as the recently proposed nonperturbative renormaliza¬ 
tion group approach of Refs. [12, HH. Our result further 
gives the full spacetime structure of the correlator for 
noncoincident points at large spacetime separations. To 
the best of our knowledge, this has never been computed 
before. It would be of definite interest to perform such 
a calculation in, say, the stochastic effective theory for 
a comparison. More generally, it would be interesting to 


make a precise link between the present approach and the 
stochastic theory, e.g., along the lines of Refs. [13, [H, [111 . 

Interesting extensions of the present work include the 
generalization of the techniques developed here to other 
types of (e.g., fermionic or gauge) fields and/or to theo¬ 
ries with more complicated (e.g., derivative) interactions. 
Also of interest is the study of other spacetime geome¬ 
tries, such as quasi-de Sitter space and general cosmo¬ 
logical spacetimes, or black hole geometries, of interest 
for the question of radiative corrections to the Hawking 
radiation (s^ . 


Appendix A: Self-energy at NLO 

Here, we compute the nonlocal contribution (1171) to 
the NLO self-energy. As explained before, the multiloop 
nature of the NLO approximation can be captured by 
the function I{p,p'). The latter resums the infinite se¬ 
ries of bubble diagrams and eventually takes a similar 
form as a propagator with modified infrared exponent; 
see Eqs. and (l25l) [we refer the reader to Ref. [Hj 
for details concerning the calculation of the functions 
If{PtP') and Ip{p,p')\. The remaining momentum (loop) 
integral in Eq. uni), can thus be treated using the methods 
developed in Ref. [gO]- For external momenta < p, 
it is dominated by internal (loop) momenta such that 
QP' ^ M, which imply rp, rp' < p. We can effectively 
limit the momentum integration to a sphere of radius 
min(/r/p, p/p') and use the leading infrared behaviors 
(1771) ((^ of the integrands. 

We first evaluate the statistical part of the self-energy 
Tip- In this case the integrand is proportional to 
— p®/p^/4 ~ Here, the neglect of the 

p-components is justified by the infrared power laws in 
Eqs. (1771) and (1771) . which results from the amplification 
of infrared fluctuations. The p-components are not am¬ 
plified. We thus get 


^F^(P,P') 


XFjpp f 1 

3A(pp')^/^+^ Jq g2l/j,2(l/-e) ’ 


(Al) 


where /3 = iz —27 = v — 2£ and r = |e-f q|, with e an arbi¬ 
trary unit vector. The remaining integral is rapidly con¬ 
vergent such that one can safely set the upper bound to 
infinity. It can then be evaluated, e.g., using the method 
of Feynman parameters. 


Hd r(e)r( 27 )r(zz + e)r(i7-2e) 


q2vp2{v-e) 2(27r)^ r(£-|-27) r(^)r(zz — e) 




(27r)« 


1 


I 

47 


(A2) 


This is a usual feature of the strong (classical) field regime 
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where we have neglected relative corrections of C>(e, 7 ) in 
the last line. We finally get 


^F^iP,p') 


o'pFu 

(pp/)3/2+^ ’ 


(A3) 


where ap is defined in Eq. (l29l) . 

The integrand for the spectral component Sp involves 
the combination For p < p', it reads 


, f 1 1 ) 

•^|q|<^l \P/ iPPT~‘' j 

(A4) 


The remaining integrals are trivially performed, 


i|q|<M/p' {27rY{d-2a) \p' 


(A5) 


Repeating the same calculation for p > p', we finally 
obtain 




A E^TTp fid / p' ) p' ) I 

3N ^pp>)l {27 t)A 2e 4j j 


(pp')3/2 


-p2jj 
' V 



(A 6 ) 


where we have neglected relative corrections of 0 (e, 7 ) in 
the constant factors^ ^ and where we have used V^'^{x) ~ 
Vu- 2 ri{x) = Vv- 2 e{x) Fi [x) for |x| > 1 in obtaining 
the last expression. 

Equations (|A3(1 and (IA 6 I 1 are the results quoted in the 
text, Eqs. and (| 2 ^ . 


Appendix B: Propagator in spacetime 

Here, we consider the spacetime structure of the NLO 
propagator obtained in the present work, Eqs. dsni)-®- 
The latter can be written as a linear superposition of two 
free massive propagators with masses m± <C 1. We thus 
consider first the spacetime structure of the propagator of 
a free light scalar field in its CTBD vacuum. In general, 
this can be emressed exactly in terms of a hypergeomet¬ 
ric function [ 1 , but the case of a light field yields a 
much simpler form, both in the coincidence limit and at 
large spacetime separations, that we now discuss. 


As mentioned in the text, it is justified to neglect such corrections 
in numerical factors [for instance, = 1 + 0{e)], but not in the 
momentum dependences since the functions Ei? and Sp are to be 
involved in convolution integrals involving large values of, e.g., 

ln(p/p')- 


The statistical propagator of a scalar field with mass 
m -C 1 receives contributions ~ Ijm^ which are gener¬ 
ated by the infrared power law (l2^ . Let us first consider 
the case of large spacelike separations, with spatial co¬ 
moving distance |X — X'p -b rj'^. In that case, the 
oscillating phase in Eq. ([9|) effectively cuts off the momen¬ 
tum integral at A"|X — X'| < 1, and it is fully justified 
to use the leading infrared behavior of the correlator for 
|Ar? 7 |, |Ar? 7 '| < 1, given by Eq. ([22]). Using Eqs. dH), (fTOl) . 
and ( 1 ^ . we have 


{pv'rF, 


riD+im'^ 


d^K e*K (x-x') 

( 2 ^ 

m' 


IX-X' 


(Bl) 


where v = d/2 — e, with e ~ m?/d^ and where we used 


d'^K e 


iK X 


(27r)d 


1 r(£) / 2 \ 

(47r)d/2 r(d/2 - e) V|X|y 
1 1 
(47r)d/2r(d/2) e|X|2e' 


(B2) 

(B3) 


In the case of large timelike separations |X — X'p <C 
^2 _|_ ^/ 2 ^ phase factor e®x (x-x ) ~ x the relevant 
integration region and the integral in ([9]) must effectively 
be cut off at K < min for the use of the leading 

infrared form ( 1221 ) —which yields the Xjmd' contribution 
to the correlator in spacetime—to be justified. For large 
time separation, i.e., |t — t'\ = 1 ^( 77 / 77 ')! 1 , we can 

replace the upper bound by , and we obtain 


F^{x,x') 


1 


Vtu+iw? + rj' 


7777 


1 


Vto+im? 


(B4) 

A2e 


where we used 

r d'^K 0 (A - K) Qd A2^ 1 

J {2ttY ~ ^2TrY 2e ~ {4nY/^r{d/2)~' 

(B5) 

Of course, the correlator is a function of the de Sitter 
invariant 


-772 - 77'2 -b |X - X' 


/|2 


27777' 

= — cosh(t — t') -b 


IX-X' 


J+P 


(B 6 ) 

(B7) 


The expressions (IBII) and (IB4() simply rewrite, up to rel¬ 
ative corrections 0{e), 




ix,x') 


1 

Bz)+i7772|2;|e 


for | 2 ;| » 1 


(B 8 ) 


in the relevant limits of large spacelike or timelike sep¬ 
arations. The leading l/rn^ behavior (IB 8 I) can also be 
obtained directly from the expression of the exact prop¬ 
agator in terms of hypergeometric functions . 
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Finally, we note that the expression (jB 8 l) also holds in 
the coincidence limit 2 —>■ 1 , for which 


F^{x,x)KiTf^F^ J 


f l^^l) 


1 


{2nY K'^'^ rin+im? 

(B9) 

Using the above considerations, we conclude that the 
NLO statistical propagator obtained in the present paper 
[Eq. (l30)) ] behaves as 




^D+i L 


C+ 


(BIO) 


for large spacetime separations [z] 1. In the coinci¬ 

dence limit we have 


F^^{x,x) 


1 r c+ c_ 


Ud+i ( m\ m?_ I ’ 


(Bll) 


which reproduces the second line of Eq. (1571) up to relative 
corrections 


Here, the precise numerical prefactor depends on the 
(quantum) state of the subhorizon modes. The above 
(standard) expression assumes the CTBD vacuum. 

The Langevin equation (EU can be turned into a 
Eocker-Planck equation for the probability distribution 
g{t, (f) of the stochastic process. The latter admits an 
attractor solution at late times, given by 


lim Q{t,ip) oc e 

t—¥00 


(C3) 


The late-time expectation value of any local observable 
0{ip) is then given by 

f d^u>0(w) e-^o+i^(¥’) 

^ jX‘r-L.v„i ■ («) 

In the case of interest here, where the potential is^^ 


the field variance can be written as 


(C5) 


Appendix C: The field variance 

This appendix is devoted to the calculation of the local 
field variance {ip^{x))/N at NLO in the 1/fV expansion 
in both the stochastic and the Euclidean de Sitter ap¬ 
proaches. 


1. Stochastic approach 

In the stochastic approach, the long wavelength part of 
a light quantum field is treated as a slowly evolving clas¬ 
sical stochastic field—owing to the strong gravitational 
enhancement of superhorizon modes—sourced by a white 
Gaussian random noise representing the quantum short 
wavelength (subhorizon) degrees of freedom which con¬ 
stantly cross out the horizon. Neglecting spatial gradi- 
ents, this is described by the following Langevin equation 

[HlHl 

(Cl) 

where t = — ln(— 77 ) is the cosmological time, V(ifi) is the 
relevant potential for superhorizon modes [5ll| . and the 
Gaussian random noise is characterized by the correlator 


iUtW)) = 




D+l 


-SabS{t - t'). 


(C2) 


This approach neglects the self-interactions of subhorizon modes 
as well as the interaction between sub- and superhorizon modes, 
which is similar in spirit to the approximation strategy employed 
in the present work. 


{^a^b) = --^a„lnZ(a,/3), 


(C 6 ) 


where we introduced the notations a = H£)+imjg /2 and 
/3 = f2p+iA/24 and where we defined, introducing the 
change of integration variable x = 


Z{a,P) = ^ 


(C7) 


with 


/(x) = ax + /3x^- (C 8 ) 

The 1/A^ expansion of the integral (IC7p is equivalent to 
a saddle-point expansion. The saddle point x = X is 
defined as /'(x) = 0 , which is solved as 




The LO solution is, thus, 


Xo = 


—a + +Aj3 


1 


4/3 a -I- + 4,(3 

and we shall only need the NLO correction 
1 1 


(C9) 


(CIO) 


X = Xo 


-I- 4/3 


+ 0{N-^). (Cll) 


Here, the square mass and coupling parameters are to be under¬ 
stood as the effective (renormalized) ones at the horizon scale 

1 [H. 
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The NLO result for the integral (IC7I) is given by the 
Gaussian integration around the saddle point, 

\nZ{a,P) ^ ^ /"(xo) + O (iV^) , 

(C12) 

from which it follows, using the definition f'{x) = 0, that 

The calculation of each individual part is straightforward. 
We obtain 

1^1 2 /3 \ 

N TV a 2 + 4/37 

= no+iM^ ~ N [1 + Xesr) ’ 

where we used the definitions (HiH) and (HD) in obtaining 
the second line. This exactly agrees with the result of the 
present work from the solution of the KBEs, Eqs. (l48l) 
and (ISOl) . 


2. Euclidean de Sitter space 

For light scalar fields with nonderivative interactions, 
it has been shown in Refs. [H, that the effective the¬ 
ory of the field zero mode on the Zl-dimensional unit 
sphere Sjj (Euclidian de Sitter space) is equivalent to 
the stochastic prescription (IC4I) for what concerns the 
calculation of the local field fluctuations. Still, it is in¬ 
structive to compute explicitly the variance of the zero 
mode at NLO in the 1/TV expansion directly in Euclidean 
de Sitter space. 

The line element is now given by (in Zl-dimensional 
spherical coordinates) 

ds"^ = de^+sm'^9 and. (C15) 

As in the main text, we consider the symmetric phase 
in which {ipa) = 0 and where two-point functions are 
diagonal in 0(TV) space, e.g., Gab = SabG. The Euclidean 
KBEs reads 


The NLO self-energy is given by^"^ 

E{x,x') = x)S^^^ {x, x') 

+ ^G{x, x')I{x, x'), (017) 

where the function I(x,x') resums the infinite series of 
bubble diagrams, as in Fig. [21 It solves the integral equa¬ 
tion 

I{x,x’) = U{x,x') + f n{x,y)I{y,x’), (CIS) 

with the one-loop bubble 

U{x,x') = -^G‘^ix,x'). (019) 

On the O-dimensional sphere, the field can be decom¬ 
posed as 

'fa{x) = ^ 

L 

where L = (L, L^-i, ..., Li) is a vector of integer num¬ 
bers with L > Ljj-i > ■■■ > \Li\ and where the D- 
dimensional spherical harmonics satisfy 

UoYiix) =-L{L +D-l)Yp^{x) (021) 

and are normalized as 

[ Yi{x)Y^,{x)=6j^^^,. (022) 

J X 

The zero mode is the constant Yg = 1/with 
Od+i the volume of the unit sphere Sd- 

The variance of the zero mode of a free scalar field is 
oc l/m^g. Hence, the zero mode of light fields in units of 
the sphere radius undergo strong fluctuations and must 
be treated nonperturbatively. Following Refs. [H, 1^ . 
the dominant contributions to the corresponding effec¬ 
tive theory is obtained by simply discarding all nonzero 
modes. The contributions of the latter can then be con¬ 
trolled by perturbative means. In particular, we write 

ipa{x) = (fa+ ipa(x) -)> ipa, (023) 

with ipa = Pa, 0^0 = fa:‘Po.(x)/flD+i- Accordingly, we 
only retain the constant contributions to the various two- 
point functions involved in the KBEs (|C16p . that is, 

G{x,x') = G + G{x,x') ^ G, (024) 


(-□3,-bmdg) G{x,x')+ / Y{x,y)G{y,x') = S^^\x,x'), 
Jy 

(016) 

where Da, is the appropriate Laplace-Beltrami operator, 
!y = S d^vVgiy), and d(^)(x,x') = S^^\x-x')/y/g{x). 


Here, both the local and the nonlocal contributions involve the 
full propagator G. Strictly speaking this corresponds to the 1/A^ 
exp ansion at NLO in the two-particle-irreducible formalism [79 . 
1^ . The standard NLO contributions are obtained from the 
latter by systematically expanding the propagator around its LO 
expression, as we do below. 
















13 


and similarly for 'S{x,x')^ I{x,x'), and n(a;,x'). 

The previous NLO expression in the zero-mode sector 
are then 


E = - 1-H — 


G 


N J no+i 37V 




(C25) 


with 


/ = n -1- = 


n 


1 — 


(C26) 


where 11 = —-IG^. We have, finally, 


A G 




i-h 


1 


D+l 


N 1 + ^Qd+iG^ 


(C27) 


Similarly the KBEs (IClfill in the zero-mode sector read 


using Eq. (IC27I) and the definitions a = ftD+im^g/2 and 
/3 = n£)_|_iA/24 introduced in the previous subsection. 

This is easily solved at NLO in 1/N, 





2 /3 A 

iV a2 + 4/3 J ’ 


(C30) 


with the LO solution 


Go = 


1 

a + \/ q ;2 _|_ 4/3 


(C31) 


A” O^j+iEG 


1 


which can be rewritten as 


2aG + 4;3G2 



N 1 -h 4/3G2 ) 


= 1 


(C28) 

As expected, Eqs. (IC30I) and (IC31|) reproduce the 
stochastic result of the previous subsection and, there¬ 
fore, the solution of the KBEs in Lorentzian de Sitter 
(C29) spacetime obtained in the present paper. 
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